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Abstract 

We consider the problem of designing an optimal quantum detector to minimize the probabil- 
ity of a detection error when distinguishing between a collection of quantum states, represented 
by a set of density operators. We show that the design of the optimal detector can be formulated 
as a semidefinite programming problem. Based on this formulation, we derive a set of necessary 
and sufficient conditions for an optimal quantum measurement. We then show that the optimal 
measurement can be found by solving a standard (convex) semidefinite program followed by 
the solution of a set of linear equations or, at worst, a standard linear programming problem. 
By exploiting the many well-known algorithms for solving semidefinite programs, which are 
guaranteed to converge to the global optimum, the optimal measurement can be computed very 
efficiently in polynomial time. 

Using the semidefinite programming formulation, we also show that the rank of each optimal 
measurement operator is no larger than the rank of the corresponding density operator. In 
particular, if the quantum state ensemble is a pure-state ensemble consisting of (not necessarily 
independent) rank-one density operators, then we show that the optimal measurement is a 
pure-state measurement consisting of rank-one measurement operators. 



Index Terms — Quantum detection, semidefinite programming, duality. 
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1 Introduction 

In a quantum detection problem a transmitter conveys classical information to a receiver using a 
quantum-mechanical channel. Each message is represented by preparing the quantum channel in 
a quantum state represented by a density operator, drawn from a collection of known states. At 
the receiver, the information is detected by subjecting the channel to a quantum measurement in 
order to determine the state prepared. If the quantum states are mutually orthogonal, then the 
state can be determined correctly with probability one by performing an optimal orthogonal (von 
Neumann) measurement [ffl]. However, if the given states are not orthogonal, then no measurement 
will distinguish perfectly between them. Our problem is therefore to construct a measurement that 
minimizes the probability of a detection error. 

We consider a quantum state ensemble consisting of m density operators {pi, 1 < i < m} on an 
n-dimensional complex Hilbert space 7i, with prior probabilities {pi > 0, 1 < i < m}. A density 
operator p is a positive semidefinite (PSD) Hermitian operator with Tr(p) = 1; we write p > 
to indicate p is PSD. A pure-state ensemble is one in which each density operator pi is a rank- 
one projector |0j)(<^j|, where the vectors \4>i), though evidently normalized to unit length, are not 
necessarily orthogonal. 

For our measurement we consider general positive operator- valued measures |2|, ||] , consisting of 
m PSD Hermitian operators {Hj, 1 < i < m} that form a resolution of the identity on TL. A pure- 
state measurement is one in which each measurement operator is a rank-one operator^ \[ii)(pi\, 
where the vectors \pi) are not necessarily orthogonal or normalized. An orthogonal measurement 
(i.e., a von Neumann measurement) is one in which the measurement operators Hi are mutually 

1 In this paper, when we say rank-one operator we mean an operator that can be expressed in the form Hi — |^tj)(/x, | 
for some £ TL. Note, however, that may be equal to in which case the operator actually has rank zero. 
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orthogonal projection operators. 

Necessary and sufficient conditions for an optimum measurement minimizing the probability of 
a detection error have been derived [|], S|. However, except in some particular cases Q |6|, ^, ||] , 
obtaining a closed-form analytical expression for the optimal measurement directly from these 
conditions is a difficult and unsolved problem. Thus in practice, iterative procedures ju]] or ad-hoc 
suboptimal measurements are used. A detection measurement that has many desirable properties 
and has been employed in many settings is the least-squares measurement H, also known as the 
square-root measurement [jn], Q . 

Holevo § derives the necessary and sufficient conditions by considering infinitesimal transforma- 
tions of the measurement operators Ilj that preserve their character as elements of a measurement. 
The drawback of this approach is that it does not readily lend itself to efficient computational 
algorithms. Yuen et al. || use the principle of duality in vector space optimization to derive the 
same necessary and sufficient conditions. Specifically, they show that the problem of finding the 
measurement that minimizes the probability of a detection error can be formulated as a generalized 
linear programming problem, with the positive orthant being replaced by the positive cone of PSD 
matrices. Although their approach leads to the same conditions derived by Holevo their ap- 
parent suggestion that this formulation produces a standard finite-dimensional linear programming 
problem is not correct, because the cone of PSD matrices cannot be described by a finite set of 
linear inequalities. 

In this paper, we derive the necessary and sufficient conditions for an optimal quantum measure- 
ment in a self-contained manner, again by exploiting duality arguments. The primary advantage 
of our formulation is that it readily lends itself to efficient computational methods. Specifically, 
we show that the optimal measurement can be found by solving a standard convex semidefinite 
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program with n 2 variables, followed by the solution of a set of linear equations or, at worst, a 
standard linear programming problem. By exploiting the many well-known algorithms for solving 
semidefinite programs [13, 14, 15, 16], the optimal measurement can be computed very efficiently in 
polynomial time. Furthermore, in contrast to the iterative algorithm proposed by Helstrom |l^] for 
solving the quantum detection problem, which is only guaranteed to converge to a local optimum, 
algorithms based on semidefinite programming are guaranteed to converge to the global optimum. 

After a statement of the problem in Section [2], we derive in Section ^ the necessary and sufficient 
conditions for the optimal measurement that minimizes the probability of a detection error, by 
formulating our problem as a semidefinite program. Using this formulation, in Section |^ we prove 
that if the quantum state ensemble is a pure-state ensemble consisting of rank-one density operators 
Pi = \4'i){4 > i\j then the optimal measurement is a pure-state measurement consisting of rank-one 
measurement operators IT = \ni)(fii\. This generalizes a previous result by Kennedy jL7|], which 
establishes that for linearly independent vectors the optimal measurement is a (necessarily 
orthogonal) pure-state measurement. We also show that for a general quantum state ensemble, 
the rank of each optimal measurement operator Ilj is no larger than the rank of the corresponding 
density matrix pi. In Section [5] we consider efficient iterative algorithms that are guaranteed to 
converge to the globally optimum measurement. 

Throughout the paper we use the Dirac bra-ket notation of quantum mechanics. In this notation, 
the elements of TL are "ket" vectors, denoted, e.g., by \x) G TL. The corresponding "bra" vector (x\ 
is the conjugate transpose of \x). The inner product of two vectors is a complex number denoted 
by (x\y). An outer product of two vectors such as \x)(y\ is a rank-one matrix, which takes \z) G TL 
to (y\z)\x) G TL. 
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2 Optimal Detection of Quantum States 

Assume that a quantum channel is prepared in a quantum state drawn from a collection of given 
states. The quantum states are represented by a set of m PSD Hermitian density operators 
{Pijl < i < m } on an n-dimensional complex Hilbert space H. At the receiver, a measurement 
is constructed, comprising m PSD Hermitian measurement operators {Ilj, 1 < i < m} on The 
problem is to choose the measurement operators to minimize the probability of detection error, 
i.e., the probability of incorrect detection of the transmitted state. 

We assume without loss of generality that the eigenvectors of the density operators {pi, 1 < i < 
m} spanf] 7i. In this case, to constitute a measurement, the measurement operators IT must satisfy 

m 

1=1 

where I is the identity operator on TC. 

We seek the PSD measurement operators {IT, 1 < i < m} satisfying (|l|) that minimize the 
probability of detection error, or equivalently, maximize the probability of correct detection. Given 
that the transmitted state is pj, the probability of correctly detecting the state using measurement 
operators {IT, 1 < i < m} is Tr(pjHj). Therefore, the probability of correct detection is given by 

m 

P d = Y,PiMpitti), (2) 

1=1 

where pi > is the prior probability of pi, with Y^iPi = Denoting by B the set of Hermitian 

2 Otherwise we can transform the problem to a problem equivalent to the one considered in this paper by refor- 
mulating the problem on the subspace spanned by the eigenvectors of {pi, 1 < i < m}. 
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operators on TL and denning p\ = piPi, our problem reduces to the maximization problem 



max > Tr(/;-II,-). (3) 



subject to the constraints 



H > 0, 1 < i < m; (4) 



Denoting by A the set of all ordered sets II = {Iljl^l^IIj S B, satisfying (|j) and (|5|), and defining 
J(I1) = YliLi Tr(Pillj), we can express our maximization problem as 



maxJ(n). (6) 

neA y ' w 



We refer to A as the feasible set, and to any II £ A as a feasible point. Since A is a compact set 
and J(II) is a continuous linear functional, there exist an optimal II G A and an optimal value J 
defined by 

J = J(fi) > J(n), vneA. (7) 

Equipped with the standard operations of addition and multiplication by real numbers, B is 
an n 2 -dimensional real vector space. By choosing an appropriate basis for B, the problem of 
(||) can be put in the form of a standard semidefinite programming problem, which is a convex 
optimization problem; for a detailed treatment of semidefinite programming problems see, e.g., 



15, 16, 13]. Rather than relying on results that are scattered throughout the literature in 
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various forms, in what follows we present a self-contained and direct derivation of the necessary 
and sufficient conditions for the optimal measurement. As we will see, this derivation also leads to 
efficient methods for computing the optimal measurement in cases in which an analytical solution 
is not known. 

In the next section we derive the necessary and sufficient conditions on the measurement oper- 
ators by formulating a dual problem. The dual problem will also be used in Section || to develop 
efficient computational algorithms. 

3 Dual Problem Formulation 

Our objective is to formulate a dual problem whose optimal value serves as a certificate for J. 
Specifically, we will formulate a minimization problem of the form minx T{X) for some linear 
functional T such that for all feasible values of X € B, i.e., values of X G B that satisfy a certain 
set of constraints, and for any II 6 A, we shall have T(X) > J(II). The dual problem therefore 
provides an upper bound on the optimal value of the original (primal) problem. In addition we 
would like the minimal value of T, denoted T, to be equal to J. The equality J = T will then 
lead to conditions of optimality on the measurement operators. Furthermore, in this case, instead 
of solving the primal problem, we can find J and the optimal measurement by solving the dual 
problem, which turns out to have far fewer decision variables. 

3.1 Constructing The Dual Problem 

A general method for deriving a dual problem is to invoke the separating hyperplane theorem 
which states that two disjoint convex sets^ can always be separated by a hyperplane. We will take 
3 A set C is convex if for any x,y G C, ax + (1 — a)y £ C for all a £ [0, 1]. 
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one convex set to be the point 0, and then carefully construct another convex set that does not 
contain 0. This set will capture the equality constraints in the primal problem and the fact that 
for any primal feasible point, the value of the primal function is no larger than the optimal value. 
The dual variables will then emerge from the parameters of the separating hyperplane. 

In our problem we have one equality constraint, YliLi H = I, and we know that J > J (II). Our 
constructed convex set will accordingly consist of matrices of the form — / + YliLi H where IT G B 
and Hi > 0, and scalars of the form r—J(H) where r > J. We thus consider the (n 2 + l)-dimensional 
real vector space 

C = B xTZ = {(S,x) : S G B, x G K}, 
where 1Z denotes the reals, with inner product defined by 



((W,y),(S,x))=Tv(WS)+yx. (8) 

Note that since W, S G B, Tv(WS) G K. 
We now define the subset Q of C by 

n= j(- / + |i n *' r -E Tr ( n ^)) : H G Ilj > 0, r ElZ, r > . (9) 

It is easily verified that Q is convex, and g" $7. Therefore, by the separating hyperplane theorem, 
there exists a nonzero vector (Z, a) G C such that ((Z, a), (Q, b)} > for all (Q, b) G ft, i.e., 

Tr (^z (^-i + n ^ + « ^ - E Tr ( n */4)^ > o (10) 



8 



for all Ilj G B and r € 1Z such that IL; > 0, r > J. It will turn out that the hyperplane parameters 
(Z, a) define the optimal dual point. We first show that these parameters have to satisfy certain 
constraints, which lead to the formulation of the dual problem. 
Note that @ with EE* = 0, r -> J implies 

aJ>Tr(Z). (11) 

Similarly, ( |10D with r = J + 1, Hj = for j ^ i, = txx' where x G 7?. n is fixed and t — » +cxd 
yields a;'(Z — ap'^x > 0. Since x and i are arbitrary, this implies 

Z>api, l<i<m. (12) 



With Hi = 0, r -> +oo, (gg) implies a > 0. If a = 0, then @ yields Z > 0, and (y) yields 
> Tr(Z), which together means Z = 0. However, this would contradict the assumption that 
(Z, a) / 0. Therefore we conclude that a > 0, and define X = Z/a. Then (|IT|) implies that 



TpT) < J, (13) 



where T(X) = Tr(X), and (12) implies that X > p[ for 1 < i < m. Let V be the set of X £ B 



satisfying X > pL 1 < i < m. Then for any X G T, II E A, we have 



r(x) - J(n) = J] tv (n,(x - ^)) > o. (14) 



i=l 



Since I6T, from (||) and @ we conclude that T(X) = J. 
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Thus we have proven that the dual problem associated with (||)-(|5|) is 



minT(X), 
xeB 



(15) 



where T(X) = Tr(X), subject to 



X > p'i, l<i<m. 



(16) 



Furthermore, we have shown that there exists an optimal X € T and an optimal value T defined 
by 

f = T(x)<T(x), vier, (17) 



such that 



T = J. 



(18) 



3.2 Optimality Conditions 

Let Hi denote the optimal measurement operators that maximize @ subject to @) and <^j, and 
let X denote the optimal X that minimizes (15) subject to (0). Then from (18) it follows that 



(19) 



t=l 



Since X > p\ and IT > 0, (^9|) is satisfied if and only if 



[X - p'i)% = Ui(X -p'i) = 0, 1 < i < m. 



(20) 



Once we find the optimal X that minimizes the dual problem (|T^), the constraint (2C) is a 
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necessary and sufficient condition on the optimal measurement operators Ilj. We have already seen 
that this condition is necessary. To show that it is sufficient, we note that if a set of measurement 
operators 11* satisfies @, then XX 1 1* (^(X - pf^j = so that J (IT) = T(X) = J. 

Since the dual problem involves many fewer decision variables than the primal maximization 
problem, it is advantageous to solve the dual problem and then use ( |20| ) to determine the optimal 
measurement operators, rather than solving the primal problem directly. In Section |H| we develop 
efficient algorithms that follow this strategy. 

Using @, (p0|) and leads to the conditions 



Y,p[% = Y^Tkfa (21) 

1=1 1=1 

m 

^f/fii > Pp l<i<m. (22) 



i=l 



Thus, any optimal measurement II = {Hj}^ must satisfy (^) and (p2|). These conditions are also 
derived in j|, ||]. However, as noted in the Introduction, the approach taken here lends itself to 
fast iterative algorithms, as we will see in Section [El, and also provides additional insight into the 
optimal measurement operators, as we show in Section |j. 

In || it was established that the conditions ( |2~I| ) and (22) together with (Q) and (|5|) are also 
sufficient. For completeness, we repeat the argument here. Suppose that the measurement operators 
Ui satisfy © and (||). Then X = Yd=i ^iPi G r - lt then follows from @ that for any set of 



measurement operators Ilj £ A, 



J2 Tr (n iP < Tr(X) = £ Tr(n^) (23) 



i=l i=l 
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with equality for H = Ilj. Therefore the measurement operators Ilj are optimal. 
We summarize our results in the following theorem: 

Theorem 1. Let {pi, 1 < i < m} denote a set of density operators with prior probabilities {pi > 
0, 1 < i < m}, and let {p[ = PiPi, 1 < i < m}. Zei A denote the set of all ordered sets of Hermitian 
measurement operators H = {Hj}^ that satisfy Hi > and YliLi H» = I> an d ^ ^ denote the set 
of Hermitian matrices X such that X > pi, 1 < i < m. Consider the problem maxrieA J(R) an d 
the dual problem minxgr^(^); where J(H) = YlTLi Tr(p'JIi) and T(X) = Tr(X). Then 

1. For any X £ T and U G A, T(X) > J{Ii). 

2. There is an optimal H, denoted II, such that J = J(H) > J(n) for any H S A; 

3. There is an optimal X , denoted X, such that T = T(X) < T(X) for any X £ T; 

4. f = J; 

5. Given X, a necessary and sufficient condition on the optimal measurement operators Hi is 
{X - p^fii = 0, 1 < % < m. 

4 Rank- One Ensembles 

Suppose now that the density operators pi are rank-one operators of the form pi = \<pi)((j)i\ for some 
\<pi) G 7i. In this case, it seems intuitively plausible that the optimal measurement will consist of 
rank-one measurement operators of the form Ilj = \pi)(pi\ for some \pi) G 7i. 

There are some particular cases in which an analytical solution to the quantum detection prob- 
lem is known ||, ||, 0, ^, || . In all of these cases, when the density operators are rank-one operators, 
the optimal measurement also has rank one. In the special case in which the vectors \4>i) are linearly 
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independent, Kennedy (T^j showed that the optimal measurement is always a rank-one measure- 
ment. However, this implication has not been proven in the general case. Using the conditions 
for optimality we derived in the previous section, we now prove this implication for an arbitrary 
rank-one ensemble. 

We have seen that the optimal measurement operators Ilj can be determined by solving fl20|), 
where X is the optimal matrix that minimizes (|i"5|) subject to Thus the measurement oper- 

ators Ilj must lie in the null space of X — p' if denoted Af(X — p^), and consequently rank(ILj < 
dim(^(X-p0). 

Since X > pi, 1 < i < m, it follows that X is positive definite on 7i. Indeed, since the 
eigenvectors of the matrices pi span TC, for any h £ Ti there exists an i such that > 0, 

which implies that (h\X\h) > for any h E TL, so that M(X) = {0}. Now, for any two matrices Z\ 
and Zi, rank(Zi + Z2) > rank(Zi) — rank(Z2), so that 



dim(AA(Zi + Z 2 )) < dim(AA(Zi)) + rank(Z 2 ). (24) 



With Zi = X and Z 2 = -d, <M) yields 



dim (M{X - p'M < rank(p^) = rank(pi), (25) 



and 

rank(fLj < dim (aT{X - p-)) < rank(pi), 1 < i < m. (26) 
In the special case in which the operators pi = \4>i){(fii\ have rank-one, it follows immediately 



from (26) that the optimal measurement operators also have rank-one, so that they have the form 
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H = \fJ,i){m\ for some G H. 

If in addition the vectors {|<^i), 1 < i < m} are linearly independent, then the vectors 1 < 

i < m} must also be linearly independent since Y^iLi is equal to the identity on TL, where 

now H is the m-dimensional space spanned by the vectors \4>i)- Then, for 1 < j < m, 

m 

\n) = ( 27 ) 

i=l 

Since the vectors are linearly independent, we must have that = 5%j so that the vectors 



are mutually orthonormal. We therefore recover the statement by Kennedy [17], that for a 
pure-state ensemble with linearly independent vectors, the optimal measurement is an orthogonal 
pure-state measurement. 

We summarize our results in the following theorem: 

Theorem 2. Let {pi, 1 < i < m} be a quantum state ensemble consisting of density operators pi 
with prior probabilities pi > 0. Then the optimal measurement consists of measurement operators 
{Hi, 1 < i < m} with rank(Ui) < rank(pi). In particular if {pi = \4>i){4>i\,l < i < m} is a pure- 
state quantum ensemble, then the optimal measurement is a pure-state measurement consisting of 
measurement operators of the form {Hi = \m)(m\, 1 < i < m}. 



5 Computational Aspects 

In the general case there is no closed-form analytical solution to the maximization problem (||) 
or the minimization problem (|i~5|). However, since ( |l5|) is a convex optimization problem, there 
are very efficient methods for solving (|l^). In particular, the optimal matrix X minimizing Tr(X) 
subject to ( p!B[ ) can be computed in Matlab using the linear matrix inequality (LMI) Toolbox. A 
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convenient interface for using the LMI toolbox is the Matlab package^] IQC/3. 

Once we determine X, the optimal measurement operators IT can be computed using (EO), (0) 
and (||). Specifically, from (20) and (0) it follows that Ilj can be expressed as 



Ui = ^2a ij \q ij )(q ij \, (28) 
3=1 

where aij > 0, t < rank(/?j), and the vectors \qij), 1 < j < t are orthonormal and span J\f(X — /?•)• 
To determine the vectors \qij) we may use the eigendecomposition of X — p\. 
To satisfy @ we must have 

5^5^0tjlftj>(«jl = 7 - ( 29 ) 

i=l j=l 

Let |e) = vec(J) and \yij) = vec(\qij) (qij\), where \v) = vec(V) denotes the vector obtained by 



stacking the columns of V. Then we can express (29) as 



Y\a) = |e), (30) 

where Y is the matrix of columns \yij) and \a) is the vector with components a^. If the matrix Y 
has full column rank, then the unique solution to (|30| ) is 



\a) = (Y*Y)- L Y*\e). (31) 

In the general case, Y will not have full column rank and there will be many solutions \a) to 
(|30|). Each such vector defines a corresponding set of optimal measurement operators IT via (|28|) . 



4 This software was created by A. Megretski, C-Y. Kao, U. Jonsson and A. Rantzer and is available at 
http : //www .mit . edu/ cykao/home .html. 
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To find a unique solution we may seek the vector^] \a) > that satisfies (|30|), and such that 

=1 dij is minimized. Our problem therefore reduces to 



min (l\a) 



(32) 



where |1) denotes the vector with components that are all equal 1, subject to 



Y\a) = |e); 
la) > 0. 



(33) 



The problem of (p2|)-(P^) is just a standard linear programming problem that can be solved very 



efficiently using standard linear programming tools [19], for example the LMI toolbox in Matlab. 
5.1 Example 

We now consider an example illustrating the computational steps involved in computing the optimal 
measurement. 

Consider the case in which the ensemble consists of 3 rank-one density operators pi = \<f>i)((j>i\, 1 < 
i < 3 where 

1 i 1 

(34) 





1 




1 







1> = 





■ '^7! 


1 


, \<h) = 


1 







with prior probabilities 



pi = 0.1, p 2 = 0.6, p 3 = 0.3. 



(35) 



3 The inequality is to be understood as a component-wise inequality. 
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To find the optimal measurement operators, we first find the optimal matrix X that minimizes 
Tr(X) subject to X > p\ with p\ = PiP\. The matrix X is computed using the IQC toolbox on 
Matlab. To this end we generate the following code: 



» abst_init_lmi 

» X=symmetric (2) ; 

» X>pl*Rl; 

» X > p2 * R2; 

» X > p3 * R3; 

» lmi_mincx_tbx (trace (X) ) ; 

» X=value(X) 

The optimal X is given by 



% Initializing the LMI toolbox 

% Defining a symmetric 2x2 variable X 

% Imposing the inequality constraints: 

% Here pi = pi,p2 = p2, p3 = ps and 

% Rl = pi,R2 = p 2 ,R3 = p 3 

% Minimizing Tr(X) subject to the constraints 

% Getting the optimal value of X 



X 



0.352 0.217 
0.217 0.434 



(36) 



Using the eigendecomposition of X — p\ we conclude that, as we expect from Theorem ^, Af(X — p[) 
has dimension 1 for each i and is spanned by the vector \qi) where 



\Qi) 



-0.833 
0.554 



\12/ 



1 

V2 



0.850 




-0.525 




> l<?3> = 




0.527 




0.851 



(37) 



The optimal measurement operators are therefore given by II j = a,i\qi)(qi\ = \pi){p,i\ with \pi) 
J7ii\qi} and a* denoting the ith component of \a). From p0) \a) must satisfy 



0.693 


0.722 


0.276 


-0.461 


0.448 


-0.447 


-0.461 


0.448 


-0.447 


0.306 


0.278 


0.724 









1 




a\ 


















a3 




1 



(38) 
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Since the matrix in ( |38| ) has full column rank, there is a unique solution 



0.007 
0.999 
0.994 



The optimal measurement vectors are then given by = ^/(ii\qi) which yields 



-0.067 




0.849 




-0.524 




» |A*2> = 




» |/*3> = 




0.046 




0.527 




0.849 



(39) 



(40) 



We can immediately verify that the measurement operators H = with given by 

( (40| ) together with X given by (^6|) satisfy the necessary and sufficient conditions (||) , (|5|) and pC 
Furthermore, we have that the probability of correct detection is given by 



Tr(X) = f^PiM^Pi) = 0-78. 



(41) 



i=l 



In Fig. |T] we plot the weighted state vectors |^) = ^il^i) given by ([34]) and (|35|), together with 
the optimal measurement vectors given by fl40|). For comparison, we also plot the least-squares 
measurement vectors \xi) which are given by || 



\ X i) = (**T 1/2 |Vi 



(42) 



where \& is the matrix of columns \ifji) and (-) 1//2 is the unique symmetric square root of the 
corresponding matrix. Note, that since the vectors span Tl, \I"I/* is invertible. The probability 



18 



of correct detection using the least-squares measurement vectors is ^Ci^iP«KXi|^i)| 2 = 0.71. As 
we expect, this probability is smaller than the probability of correct detection using the optimal 
measurement vectors which from (41) is equal to 0.78. 



|M3> |X3) 




Figure 1: Illustration of the optimal measurement vectors. The weighted state vectors are \ipi) = 
■s/p~i\4>i) where the vectors \4>i) and the probabilities pi are given by (|3^ ) and (|35l) , respectively. The 
optimal measurement vectors \ni) are given by (|40|). The least-squares measurement vectors \xi) 
are plotted in dashed lines for comparison, and are given by (42). 
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